arXiv: 1507.07562vl [astro-ph.SR] 27 Jul 2015 


Astronomy & Astrophysics manuscript no. bragaglia'rv © ESO 2015 

July 29, 2015 


NGC 6139: a normal massive globular cluster or a first-generation 
dominated cluster? Clues from the light elements * ** 

A. Bragaglia 1 , E. Carretta 1 , A. Sollima 1 P. Donati 1 - 2 , V. D’Orazi 3,4,5 , R.G. Gratton 3 , S. Lucatello 3 , and C. Sneden 6 

1 INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1,1-40127 Bologna, Italy 

2 Dipartimento di Fisica e Astronomia, Universita di Bologna, viale Berti Pichat 6,1-40127 Bologna, Italy 
1 INAF-Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5,1-35122 Padova, Italy 

4 Monash Centre for Astrophysics, School of Physics and Astronomy, Monash University, Melbourne, VIC 3800, Australia 

5 Department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109, Australia 

6 Department of Astronomy and McDonald Observatory, The University of Texas, Austin, TX 78712, USA 


ABSTRACT 

Information on globular clusters (GC) formation mechanisms can be gathered by studying the chemical signature of the multiple populations 
that compose these stellar systems. In particular, we are investigating the anticorrelations among O, Na, Al, and Mg to explore the influence 
of cluster mass and environment on GCs in the Milky Way and in extragalactic systems. We present here the results obtained on NGC 6139 
which, on the basis of its horizontal branch morphology, had been proposed to be dominated by first-generation stars. In our extensive study 
based on high resolution spectroscopy, the first for this cluster, we found a metallicity of [Fe/H]=-1.579 ± 0.015 ± 0.058 (rms=0.040 dex, 45 
bona fide member stars) on the UVES scale defined by our group. The stars in NGC 6139 show a chemical pattern normal for GCs, with a 
rather extended Na-0 (and Mg-Al) anticorrelation. NGC 6139 behaves like expected from its mass and contains a large fraction (about two 
thirds) of second-generation stars. 

Key words. Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: globular clusters - Galaxy: globular clusters: individual: 
NGC 6139 


1. Introduction 


Once considered as a good example of simple stellar popula¬ 
tions, Galactic globular clusters (GCs) are currently thought 
to have formed in a complex chain of events, which left a 
fossil record in their chemical composition (see e.g., the re¬ 


view by Gratton et al., 2012). Photometrically, GCs often ex¬ 


hibit spread, split and even multiple sequences that can be 
explained by different chemical composition among cluster 
stars, in parti c ular of l ight elements li k e He, C, N, O (e.g. 


Carretta et all 20 lbl Sbordone et al. . 2011 : Milone et al 


2012al : IPiotto et al. . 120151) . Our FLAMES survey of more than 


20 Milky Way (MW) GCs (see ICarretta et all f2009al l2014hl 
and references therein) combined with literature data, demon¬ 
strated that most, p erhaps all, GCs host multiple stellar popu¬ 
lations (see Carretta et al. . 2010a i). Variations of Na and O (and 


Send offprint requests to: A. Bragaglia, angela.bragaglia@oabo.inaf.it 
* Based on observations collected at ESO telescopes under pro¬ 
gramme 093.B-0583 and on data obtained from the ESO Science 
Archive Facility under request number 94403 

** Tables xx, yy, and zz are only available in 

electronic form at the CDS via anonymous ftp 

to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/???/??? 


sometimes of Mg and Al) abundances trace these different sub¬ 
populations. 

Our large and homogeneous database allowed us a quan¬ 
titative study of the Na-O anticorrelation. In all the analyzed 
GCs we found about one third of stars of primordial compo¬ 
sition, similar to that of field stars of similar metallicity (low 
Na, high O), i.e., belonging to the first generation (FG). The 
other two thirds have a modified composition (increased Na, 
depleted O) and belong to the seco nd generation (SG) of stars, 
polluted by the FG (see, e.g., Carretta et al. . 2009al bl). Which 
were the more massive stars of FG that produced the gas of 
modified comp o sition i s still an unsettled q ues tion, see e.g.. 


Decressin et al.l (120071) : IVentura et alJ (120011) . or iBastian et al 
(2013) for an alternative view. 


We also found that the extension of the Na-O anticorrela¬ 
tion tends to be larger for higher mass GCs and that, apparently, 
there is an observed mini mum c luster mass for appearance of 
the Na-O anticorrelation ( ICarretta et al„ 2010a[) . It is important 
to understand if this limit is due to the small statistics (less 
low-mass clusters have been studied, and only a few stars were 
observed in each). This in an important constraint for cluster 
formation mechanisms, because it indicates the mass at which 
we expect that a cluster is able to retain part of the ejecta of the 
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Fig. 1 . (a) CMD of NGC 6139, with its quite short blue HB. In grey all stars, in red all targets observed with FLAMES. The 
seven UVES stars are indicated by blue open squares; the stars with RV from Saviane et al. ( 20121) by blue crosses, (b) The non 
member stars (and the star on the HB) are indicated by red open squares, the stars candidate members on the basis of their RV by 
filled blue squares. The member stars define a reasonably tight RGB and we indicate with green circles the probable AGB stars 
(see Sec. 3.1) 


FG, hence to show the Na-O signature (the masses of the orig¬ 
inal clusters are expected to be much higher than the present 
ones, since the SG has to be formed by the ejecta of the FG). 

Variations in Na, O, and He are tied, but do not tell us ex¬ 
actly the same story. The effects of increased He are visible 
in colour-magnitu de diagrams (CMD) in the main sequence 


the horizontal branch phase (HB 

e.g., 

D’Antona et al. 2005; 

Gratton et al. 201(1 Milone et al. 

2014 

). Based on the pos- 


sibility of reproducing t heir HBs with a single He value. 


(e.g., Prantzos et all 20071 for an application to the GC pe¬ 
culiar chemistry), i.e. by more massive stars, should follow 
the He enrichment better than Na (and O). We can then ex¬ 
pect to find Na and O variations without significant He en¬ 
hancement; this would produce short HBs (and no variations in 
Al){Jus is what seems to happen, for instance, f or NGC 6397 
and NGC 6838 ( Carretta et all 2009bl) . for which Gratton et al 


Caloi & D’Antonal ( 201 I f. proposed that some GCs are com¬ 


posed by a single generation (FG-only) or predominantly by 
a single generation (mainly-FG). It would be important to 
measure their Na, O, and Al abundances to clarify the is¬ 
sue. In fact Al, produced only at higher temperatures than Na 


( 2010h found n o necessity o f He di spersion to explain their 
HBs. Hovever, Milone et al. (2012b[), studying the main se¬ 
quence of NGC 6397, found an internal variation in Y of about 
0.01, small, but exceeding their measurement error. In sum¬ 
mary, we would be in presence of FG and SG stars even with¬ 
out significant He variations; this is important to understand the 
cluster formation and early intracluster gas pollution. 
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After studying the high-mass clusters, we begun a system¬ 
atic study of low-mass GCs, FG-only/mainly GCs, and high- 
mass and old open clusters (OCs). Our goal is to empirically 
find the mass limit for the appearance of the Na-O anticorre¬ 
lation and to ascertain whether there are differences between 
high-mass and low-mass clusters’ properties, e.g. in the rela¬ 
tive fraction of FG and SG stars. 

We also included GCs belonging to the Sagittarius dwarf 
spheroidal (Sgr dSph) to study if there are differences in GCs 
formed in different environments (the MW and dwarf galaxies). 
In fact, GCs born in a dSph are expected to have experienced 
a milder tidal field thus possibly retaining a larger fraction of 
their original mass. Only a few old GCs in Fornax and LMC 
have their abundances derived using high-res spec t rosco py. 


These GCs also seem to host two population s dLetarte et al 


20061: Johnson et al. . 20061: Mucciarelli et al. . 20091) but the 
fractions of FG and SG stars in Fornax and LMC GCs of 
similar mass seem different. Is this again a problem of low 
statistics or is the galactic environment (a dwarf spheroidal vs 
a dwarf irregular) influencing the GC formation mechanism? 
This is outside the main theme of the present pape r and we re- 
fer the interested read er to the papers by lLarsen et al.l (120121) : 
D’Antona et alJ ( 2013 ) for a more detailed discussion. 

We have already gathered high-resolution spectra of a low- 
mass GC (NGC 6535, which will be analysed in a forthcom¬ 
ing paper), two high-mass OCs (Berkeley 39, NGC 6791, 


see 


Bragaglia et al. 20121 20141) . a proposed FG-only GC 


(NGC 6139, this paper), and Terzan 8 (ICarretta et al., 2014aj), 


Tautvaisiene etalj |2004^ Sbordone et alJl2007 ; Cohenll2004 ) or 
to Rup 106 ( Villanova et all 2013 ). However, in Terzan 8 the 
SG seems to be a minority component, contrary to what occurs 
for high-mass GC s. The Na-O anticorrelation had never been 
observed in OCs ( de Silva et al. . 20091 MacLean et al. . 2014 ) 
and we confirmed thi s usi ng large samples of stars both for 
Berkeley 39 ( Bragaglia et al .. 20 121) u sing FLAMES spectra 
and NGC 6791 dBragaglia et all E()l4 using Hydra@WIYN 
and HIRES @Keck spectra. In particular, for NGC 6791, where 
Geisler et al. ( 2012 ) claimed to have found some variations in 
Na, O, we did not detect any evidence of this trend. Our find 


ings a r e corroborated by indep endent works, see ICunha et al 


lar, it has never been studied with high-resolution spectroscopy 
before. 

Hazenl ( 1991 ) found ten variable stars in NGC 6139, five 


of which are RR Lyrae stars. The four RRab seem to indi¬ 
cate a Oosterhoff type II. A colou r-ma gnitude diagram (CMD) 
was presented by Samus’ et al. ( 1996 ): they employed photo¬ 
graphic B and V plates and their photometry reaches only the 
red giant branch (RGB) and HB. They found the cluster quite 
metal-poor (close to [Fe/ H1=-2) from the slope of the RGB, 


belonging to the Sgr dSph. In Terzan 8 we see some indication 
of a SG, at variance with other low-mass Sgr GCs (Ter 7, Pal 12 


and highly reddened. Zinn & Barnes! (19981) used VI photome¬ 
try barely reaching the main sequence turn-off, noted the high 
and differential reddening (and corrected for a gradient), and 
determined [Fe/H |=— 1.71 ± 0.20 and a mean reddening of 
E(V-I) - 1.03 ±0.04, correspon ding to E(B-V) - 0.7 6+0.03. 
Similar results were reached by Ortolani et a l. ( 19991) on the 
basis of VI photometry and Davidge ( 19981) who used near- 
infrared data and who also noted that the differential reddening 
is not a significant problem, since the RGB sequence becomes 
very well defined once the field stars contamination is statisti¬ 
cally removed. Finally, Piotto et all ( 2002h presented WFPC2 
photometry, as part of their HST snapshot programme (74 MW 
GCs observed with the F439W and F555W bands) the RGB 
and HB are very well defined, and the main sequence turn-off 
is better defined than in previous works. Their data are pub¬ 
licly available and we used them for the present work (see next 
Section) 

The only determi natio n of metallicity based on spec¬ 
troscopy is by Saviane et al. (2014), who used spectra at reso¬ 
lution R ~ 2500 in the region of the infrared Ca n triplet (CaT). 
NGC 6139 is one of the 20 GCs in the paper. They observed 
19 stars, 15 of which were considered members (we used 
this information to select our targets, see next Section). The 
metallicity they obtained from the CaT lines is | Fe/Hl=-1.63, 
rms=0 .13 (on the metallicity scale defined in ICarretta et al 
2Q09cj). 


Caloi & D’ Antonal ( 2011 ) include NGC 6139 among their 


(2015); Boesgaard et al. (120151 ). 

We concentrate here on the possibly FG-only cluster 
NGC 6139. In Section 2 we present literature information on 
the cluster, in Section 3 we describe the photometric data, the 
spectroscopic observations, and the derivation of atmospheric 
parameters. The abundance analysis is presented in Section 4 
and a discussion on the light-element abundances is given in 
Section 5. 


2. NGC 6139 in literature 

NGC 6139 is a massive GC (its absol u te visu al magnitude, a 
proxy for mass, is My = -8.36, Harris 1996 ) located toward 
the centre of the MW, at l - 342.37°, b = 6.94°. The cluster 
has received relatively scarce attention in the past; in particu¬ 


list of candidate FG clusters. They tried to identify those clus¬ 
ters whose HB could be reproduced by a single mass (in the 
framework where the dispersion in mass of HB stars is due - 
at least in part- to variations in He content, in turn indication 
of multiple generations). A firs t clue is the_short HB (similar 
to the HB of NGC 6397, see Piotto et al. 2002 , which, how¬ 
ever has both a “normal” Na -O anticorrelation, and a small 
AY, see Carretta et al. 2009allb : Milone et al.l 2012b . respec¬ 
tively). However, given the high mass (logM = 5.58 M 0 , 
their Table 1), NGC 6139 should have produced a SG, but 
failed to lose a large fra ction of the FG stars, so i t should be 
presently FG-dominated. lCaloi & D’ Antonal (120111) do not dis¬ 
cuss in detail the case for NGC 6139, but its inclusion in the list, 
combined with the absence of previous high-resolution spec¬ 
troscopic studies made it a good target for our on-going pro¬ 
gramme. 

3. Observations and analysis 

Of the photometric data discussed in previous section, only the 
HST catalogue is publicly available. However, the very small 
field of view (FoV, about 2' side) is not well suited to the 
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Table 1. Log of FLAMES observations. 


Setup 

UT Date 
(yyyy-mm-dd) 

UT fm , 

(hh:mm:ss) 

exptime 

(s) 

airmass 

seeing 

(arcsec) 

HR11 

2014-05-10 

02:25:58.027 

3600 

1.459 

0.76 

HR11 

2014-05-10 

03:27:50.586 

3600 

1.217 

0.88 

HR11 

2014-05-10 

04:29:58.656 

3600 

1.091 

0.64 

HR11 

2014-05-10 

05:36:41.844 

3600 

1.035 

0.68 

HR11 

2014-05-10 

06:50:54.214 

3600 

1.052 

0.57 

HR11 

2014-05-10 

07:54:38.745 

3600 

1.136 

0.47 

HR13 

2014-07-22 

03:21:00.913 

3600 

1.164 

1.06 

HR13 

2014-07-22 

04:24:36.587 

3600 

1.366 

1.09 

HR13 

2014-07-25 

03:14:57.248 

3600 

1.178 

0.76 

HR13 

2014-07-26 

04:12:36.160 

3600 

1.383 

0.59 

HR21 

2014-06-03 

01:01:28.004 

3090 

1.409 

0.86 

HR21 

2014-07-03 

00:29:34.582 

3090 

1.137 

0.99 


FLAMES FoV (25' diameter). We then retrieved from the ESO 
archive some B and V filter frames acquired (under program 
68.D-0265) with the Wide Field Imager (WFI) at the 2.2m 
ESO-MPG telescope, which has a FoV of about 30' side. 

They were reduced in the standard way, correcting for 
bias and flat field using IRAF@ Stars were detected indepen¬ 
dently on the B and V frames and the instrumental magni¬ 
tudes were obtained using the Point Spread Function (PSF ) fit¬ 
ting code DAOPHOT-II/ALLSTAR (S tetson! ,! 19871119931) . We 
emp loyed the 2 Micron All Sky Survey Catalogue (2MASS, 
Skrutskie et al. 20061) and the CataXcorr codc0, developed by 


P. Montegriffo, to compute the astrometric solution and trans¬ 
form the instrumental pixel coordinates into J2000 celestial co¬ 
ordinates. The astrometric precision is about 0.2 arcsec, per¬ 
fectly compatible with the requirements of the FLAMES ob¬ 
servations. No standard stars were available, so we calibrated 
our photometry to the HST one using the stars in common. The 
final photometric catalogue will be made available through the 
CDS. 


The resulting CMD is shown in Fig. [7] As expected from its 
Galactic position, the field stars contamination is conspicuous, 
but the cluster RGB and HB are visible, especially when re¬ 
stricting to the very central region. The cluster sequences are 
also affected by differential reddening (DR), as already dis¬ 
cussed in literature. This could be relevant for the spectroscopic 
analysis, since our atmospheric parameters are derived from 
photometry. However, as already done for other difficult cases, 
resorting to optical-IR colours (in particular, V - K ), greatly 
all eviates the problem. For instance, the bulge GC GC 6441 
( Gratton et al .l 20071) has an rms scatter of 0.05 mag in E (B-V), 
resulting in a random (star-to-star) uncertainties in the effective 
temperatures of +80 K. We roughly evaluated the size of DR 
for NGC 6139 by defining a red giant branch (RGB) ridge line 
with the help of the HST photometry and projecting the can¬ 
didate members (see next section) on it along the reddening 


2 IRAF is distributed by the National Optical Astronomical 
Observatory, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 

3 http://www.bo.astro.it/~paolo/Main/CataPack.html 
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Fig. 2. Upper panel: histogram of all measured RVs (open his¬ 
togram) and of the candidate member stars only (blue filled 
histogram). Lower panel: RVs versus distance from the cluster 
centre. The candidate member stars are indicated by filled blue 
squares, the UVES targets by large red crosses, while the lines 
indicate the cluster mean RV and the ±3cr interval. 


vector. The average displacement required is 0.03 mag so we 
expect a possible error of the order of 50 K. 


3.1. FLAMES spectra 


NGC 6139 was observed _with the multi-object spectro¬ 


graph FLAMES@VLT (IPasauini et all 2002). We used the 
GIRAFFE high-resolution setups HR11, HR13, and HR21 
(R=24200, 22500, 16200, respectively) which contain two Na 
doublets, the [O i] line at 6300 A, and the Al doublet at 8772 A, 
plus several Mg lines. The observations were performed in ser- 
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vice mode; a log is presented in Table [T] The HR11 and HR13 
GIRAFFE observations were coupled with the high-resolution 
(R=47000) UVES 580nm setup (AA ^ 4800 - 6800 A) and 
the HR21 observations with the UVES 860nm setup (AA - 
6600 - 10600 A). We will use only the 580nm spectra in this 
paper, for uniformity with our previous works and because they 
have better signal to noise. 

The B , V WFI catalogue was used to select candidate tar¬ 
gets for the spectroscopic observations. As done for all other 
GCs in our sample, we chose only stars without close neigh 


hours (< 2.5")- We cross-identified the stars in Saviane et al 


( 2012 ) to allocate seven UVES fibres to high probability cluster 
members; the eigth fibre was put on an empty spot for sky sub¬ 
traction. We allocated 103 GIRAFFE fibres to RGB stars, one 
to a HB star, and 16 to sky positions. Virtually all targets are 
brighter than the RGB bump (V = 17.867 ± 0.019; lNataf et al 


RV=6.7+6.0 km s 1 in Harris ( 19961) . based on much lower 


resolution and qual ity dattfl and is in good agreement with 


Saviane et al. ( 2012f) . who give RV=34±4 km s 1 , based on 15 


20131) . The 111 stars observed are shown in Fig. |TJa) with dif¬ 


ferent symbols; their coordinates and magnitudes are given in 
Table[2](for candidate members, based on their RV) and Table[3] 
(for non members). All stars observed are within 7 arcmin form 
the centre, see Fig.[2j lower panel, i.e., well within the tidal ra¬ 
dius (10.5 armin, Harris 19961) . Almost all stars observed are 
beyond the half-mass radius (0.85 arcmin, Harris). 

The spectra were reduced using the ESO pipelines for 
UVES-FIBRE and GIRAFFE data; they take care of bias and 
flat field correction, order tracing, extraction, fibre transmis¬ 
sion, scattered light, and wavelength calibration. We then used 
IRAF routines on the 1-d, wavelength-calibrated individual 
spectra to subtract the (average) sky, shift to zero RV, and com¬ 
bine all the exposures for each star. The region near the [O i] 
line required special attention because of the strong sky emis¬ 
sion and the many absorptions. Given the low RV of the clus¬ 
ter, the HR13 exposures were scheduled when the Earth mo¬ 
tion took the sky emission farther away from the stars’ [O i] 
line. We further checked that in all these exposures the tel¬ 
luric absorptions did not fall on the stars’ line, by visual com¬ 
parison wdth_j_telfuric_template ; The RV was measured using 
DOOp ( Cantat-G audin et all 20141), an automated wrapper for 
DAOSPEC ( Stetson & Pancinol 120081) : the average, heliocen¬ 
tric value for each star is given in Tables [2] and [3 

We show in Fig. [2] the histogram of the RVs and aplotof 
RVs versus distance from cluster centre (taken by IHam si ! 19961 . 


2010 web update); the cluster signature is clear. Our is the first 
estimate of RV based on a large numebr of stars observed at 
high-resolution; we found (RV) - +28.88, rms=7.34 km s _1 
(a typical rms, for a GC). This value is compatible with the 


stars (for a direct comparison of RVs for eight stars in common, 
see Table [2]). 

We are left with 50 candidate members based on the RV, 
i.e., falling within ±3cr of the cluster average. After pruning 
the sample, the RGB of NGC 6139 is well defined and tight. 


even in presence of DR (see Fig. |TJb)). In the following we 
will discuss only them. 

Following our well tested procedure (for a lengthy descrip¬ 
tion, see e.g., Carretta et al. . 2009allb ). effective temperatures 
T e ff for our targets were derived using an average relation be¬ 
tween apparent magnitudes and fi rst-pass temperatures from 
V-K colours and the calibrations of Alonso et al . (1999if2001 ). 
This method permits to decrease the star-to-star errors in abun¬ 
dances due to uncertainties in temperatures, since magnitudes 
are less affected by uncertainties than colours. This is partic¬ 
ularly true for NGC 6139, which presents high and variable 
reddening, for which we used the apparent K magnitudes in 
our relation with 77 if because the impact of the DR on these 
magnitudes is smaller. The adopted reddening E(B -V) = 0 .75 
and distance modulus (m - M) v = 17.35 are from the Harris 
( 19961) catalogue. Gravities were obtained from apparent mag¬ 
nitudes and distancejnodulus, assuming the bolometric correc¬ 
tions from Alonso et al . (19991) . We adopted a mass of 0.85 M Q 
for all stars and Mboi,© = 4.75 as the bolometric magnitude for 
the Sun, as in our previous studies. 

We eliminated trends in the relation between abundances 


Kir 


Harrnjj 19961) refers to IWebbinkl 1 198 ill . who uses RVs from 
.ij il959l)~ derived on spectra at R ~ 200 of four stars. 


from Fe i lines and expected line strength (Magain, 1984) to 
obtain values of the microturbulent velocity v,. 

Finally, using the above values we interpolated within the 
Kuruczl ( 1993 ) grid of model atmospheres (with the option for 
overshooting on) to derive the final abundances, adopting for 
each star the model with the appropriate atmospheric param¬ 
eters and whose abundances matched those derived from Fe i 
lines. Five stars among the candidate members turned out to 
have a metallicity higher by more than 0.5 dex than the aver¬ 
age of the others; they are probably field stars and have been 
excluded from the most secure sample of RV and metallicity 
members. The atmospheric parameters for all bona fide mem¬ 
ber stars are given in Table [4] 

Six of the candidate members (two with UVES and four 
with GIRAFFE spectra) turned out too metal-poor (by about 
0.15 dex); their position in the CMD seemed to indicate that 
they could be Asymptotic Giant Branch stars (AGB) and not 
RGB. Following the same procedure used for the RGB stars, 
we derived a separate colour-temperature relation for them, 
finding temperatures higher by about 110 K. Adopting the new 
temperatures we repeated the analysis, finding a metallicity in 
very good agreement with the cluster average. These second 
values for 7' c ir and metallicity are given in Table [4] 

4. Abundances 

Beside Fe, we present here abundance of O, Na, Mg, and Al. 
The last was derived from the Al i doublet at 6696-98 A for 
stars observed with UVES, and from the doublet at 8772 A 
for stars observed with GIRAFFE. The abundance ratios for all 
elements are given in Table [3 together with number of lines 
used and rms scatter. 

The abundances were derived using equivalent widths 
(EW). We measured th e EW of iron and other elements us¬ 
ing the code ROSA ( Grattonl 1988 ) ado pting a relationship be - 
tween EW and FWHM (for details, see lBragaglia et al~ 2001 ). 
The atomic data for all the lines in the UVES spectra and in the 
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Fig. 3. Anticorrelations between light elements: Na-O (upper 
panel) and Mg-Al (lower panel). GIRAFFE stars are indicated 
by red filled circles and UVES stars by filled blue squares; up¬ 
per limits in O are shown as arrows. The errorbars are the clus¬ 
ter average internal (star-to-star) errors for each element for the 
GIRAFFE sample in red and UVES in blue (see text). In the 
upper panel we indicate the sepa ratio n between P and I, and I 
and E populations, according to Carretta et ah (2009a). 


HR11, HR13 setups and the solar reference values come from 


Gratton et al. (2003). The Na abundances were corrected for 


departure form local thermodynamical equilibrium according 


to 


Gratton et al. ( 19991) . 


The derived Fe abundances do not show any trend with T e ff. 
The average metallicity derived from stars with UVES spectra 
is [Fe/H]= -1.579 ± 0.015 ± 0.058 (rms = 0.040 dex, 7 stars) 
from neutral species, where the first error is from statistics and 
the second is systematic (see below for the computation). For 
stars with GIRAFFE spectra we found a very similar value: 
[Fe/H] = -1.596+0.006+0.042 (rms = 0.038 dex, 38 stars). The 
abundance derived from single ionized species is in very good 


agreement. We found [Fe/H]n= -1.541 ±0.016 ±0.054 (rms = 
0.043 dex, 7 stars) and [Fe/H]n= -1.579 ± 0.008 ± 0.042 (rms 
= 0.048 dex, 38 stars) for UVES and GIRAFFE, respectively. 
This supports the adopted temperature scale and gravities. 

Ours is the first high-resolution spectroscopic study of this 
cluster, so no real comparison with previous determinations is 
possible. However, the metallicity we find is i n good agree 


ment with the one based on CaT ([Fe/H |= -1.63, Saviane et al 


-1.7 


20121) and with the general finding of a low metallicity (~ 
or -2) based on photometry. 

To estimate the error budget weclosely followed the proce¬ 
dure described in lCarretta et al. ( 2009allbh . Table [6] (for UVES 
spectra) and Table [7] (for GIRAFFE spectra) provide the sensi¬ 
tivities of abundance ratios to errors in atmospheric parameters 
and EWs and the internal and systematic errors. For systematic 
errors we mean the errors that are different for the various GCs 
considered in our series and which produce scatter in relations 
involving different GCs; however, they do not affect the star-to- 
star scatter in NGC 6139. The cluster uncertainty in T e fr can be 
estimated by multiplying the slope of the relation T e ff — (V — K) 
in NGC 6139 and the uncertainty in E (V - K) (assumed to be 
0.055 mag). We also quadratically sum a contribution from a 
conservative estimate of 0.02 mag error in the zero point of 
V - K colour. The resultant uncertainty is propagated together 
with those in distance modulus and stellar mass to estimate the 
systematic uncertainty in surface gravity, while the systematic 
error in v t is simply obtained by dividing the internal error for 
the square root of the number of stars. The cluster (system¬ 
atic) error in the metallicity for NGC 6139 is then obtained 
by the quadratic sum of the above three terms (multiplied for 
the proper sensitivity) with the statistical errors of individual 
abundance determination. The sensitivities were obtained by 
repeating the abundance analysis for all stars, while changing 
one atmospheric parameter at the time, then taking the average; 
this was done separately for UVES and GIRAFFE spectra. The 
amount of change in the input parameters used in the sensitivity 
computations is given in the Table header. 


5. Light elements anti-correlations 

Our final sample consists of 45 giant stars (seven observed 
with UVES and 38 with GIRAFFE, no stars are in common). 
We could measure Na, Mg, and Al in all spectra (UVES and 
GIRAFFE), while O was measured in all seven UVES stars but 
only in 36 GIRAFFE stars (29 actual detection and seven upper 
limits). 

The resulting relations between O and Na, Mg and Al are 
shown in Fig. 0 where we clearly see the normal pattern of 
light elements abundance ratios in GCs. NGC 6139 displays 
a Na-O anticorrelation like all massive GCs studied to date. 
Interestingly, it also shows a large variation in Al coupled with 
a moderate variation i n Mg , A va riatio n in Al is not seen in all 
GCs (see Fig. 6 in Carretta et all 2009bl) ; however, also when 
Al varies the moderate variation in Mg is typical, with only 
rare e xceptions such as e.g ., NGC 2808, NGC 6752 ( CarrettaL 
2014 : Carretta et al., 2012). On the other hand, we observe a 
clear Na-Al correlation, with a Pearson correlation coefficient 
rp — 0.70, which is very significant. This confirms that the Na- 
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Table 6. Sensitivities of abundance ratios to variations in the atmospheric parameters and to errors in the equivalent widths, and 
errors in abundances for stars of NGC 6139 observed with UVES. 


Element 

Average 

T e ff 

l°gg 

[A/H] 

v t 

EWs 

Total 

Total 


n. lines 

(K) 

(dex) 

(dex) 

kms -1 

(dex) 

Internal 

Systematic 

Variation 


50 

0.20 

0.10 

0.10 




Internal 


4 

0.04 

0.04 

0.06 

0.01 



Systematic 


46 

0.06 

0.06 

0.02 




[Fe/H]i 

86 

+0.060 

-0.002 

-0.008 

-0.028 

0.011 

0.021 

0.058 

[Fe/H]n 

16 

-0.031 

+0.086 

+0.022 

-0.014 

0.026 

0.034 

0.042 

[0/Fe]i 

1 

-0.044 

+0.084 

+0.038 

+0.027 

0.104 

0.108 

0.096 

[Na/Fe]i 

2 

-0.016 

-0.038 

-0.015 

+0.022 

0.074 

0.075 

0.088 

[Mg/Fe]i 

2 

-0.016 

-0.009 

- 0.000 

+0.012 

0.074 

0.074 

0.021 

[Al/Fe]i 

2 

-0.017 

-0.008 

-0.002 

+0.026 

0.074 

0.075 

0.077 


Table 7. Sensitivities of abundance ratios to variations in the atmospheric parameters and to errors in the equivalent widths, and 
errors in abundances for stars of NGC 6139 observed with GIRAFFE. 


Element 

Average 

T e ff 

l°g£ 

[A/H] 

Vt 

EWs 

Total 

Total 


n. lines 

(K) 

(dex) 

(dex) 

kms -1 

(dex) 

Internal 

Systematic 

Variation 


50 

0.20 

0.10 

0.10 




Internal 


4 

0.04 

0.04 

0.13 

0.02 



Systematic 


46 

0.06 

0.05 

0.02 




[Fe/H]i 

35 

+0.058 

-0.001 

-0.007 

-0.025 

0.022 

0.041 

0.054 

[Fe/H]n 

2 

-0.035 

+0.088 

+0.020 

-0.007 

0.067 

0.101 

0.042 

[0/Fe]i 

1 

-0.041 

+0.081 

+0.034 

+0.023 

0.095 

0.144 

0.053 

[Na/Fe]i 

3 

-0.016 

-0.040 

-0.018 

+0.017 

0.067 

0.084 

0.043 

[Mg/Fe]i 

2 

-0.021 

-0.009 

-0.003 

+0.013 

0.095 

0.100 

0.021 

[Al/Fe]i 

2 

-0.030 

-0.006 

-0.002 

+0.019 

0.095 

0.101 

0.051 



-6 -8 -10 3.5 4 4.5 5 

M v (tot) logT™ f ax (HB) 


Fig. 4. Feft-hand panel: correlation between total absolute magnitude and the interquartile range of the [O/Na] ratio for GC in our 


_ _ _ _ _ag n _ _ 

FLAMES survey ( Carretta et al. . 2009al 201 Obi 2011 ai 20J_3b[_ 21H4bl 2015, and references there in). Right-hand panel: relation 


between the maximum temperature along the HB dRecio-Blanco et all 2006 ; Carretta et all 20071) and IQR[0/Na]. In each panel 
NGC 6139 indicated by a filled green star. The Spearman rank correlation coefficient r s and the Pearson linear correlation 
coefficient r p are listed in each panel. 
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Na and Mg-Al cycles are related, but also suggests that they do 
not occur exactly in the same pollutin g stars. 


Using the criteria defined in ICarretta et al.1 (l2009al) we 


can define the fraction of stars belonging to the primordial 
(P) and second-generation (I, E) components. The fraction of 
first-generation P stars is 26% ± 8%; the fraction of second- 
generation I stars (with moderate modification of the abun¬ 
dances) is 74% + 13%; finally, there are no second-generation 
E stars (with extremely modified composition) in our sample of 
43 stars in NGC 6139. These numbers agree with what is found 
in all massive Milky Way GCs (see e. g., Carretta et all 2010al 


2014bt Johnson & Pilachowski, 20121 and references therein); 


SG stars are presently the dominant population in them. Our 
study then doe s not confirm NGC 6139 as a mainly-FG cluster, 
as proposed by Caloi & D’Antona ( 201 ll) . 

A s it happened for the other GCs in our sample (e.g., 
Carretta et all 2009ah . where most of the abundance ratios are 


based on GIRAFFE spectra, the errors on abundances produce 
almost continuous distributions (Fig.[T|. However, a separation 
between the P and I stars is visible in the upper panel (Na and 
O), especially for the stars observed with UVES. A clearcut 
separation in groups along the Na-O anticorrelation has been 
seen only in hig h-quality and resolu tion data, such as t he UVES 
spectr a in M4 ( Marino et al. . 2008 1 and NGC 2808 ( Carretta . 
20141) . The Mg-Al distribution (Fig. [3] lower panel) looks con¬ 


tinuous, with a total excursion of about 1 dex in [Al/Fe]. 

The Na-O anticorrelation has moderate extension, with 
an interquartile range for the [O/Na] ratio ( CarrettaL 2006) 
IQR[0/Na] =0.647. This IQR makes NGC 6139 fit in the rela¬ 
tions with absolute magnitude (My, i.e., mass) and maximum 
tempe rature reached on the HB ((H B). 1R ec i o - BI a n c o et al 


20061: Carretta et al.l 20071 see Fig. [4] for a graphical represen¬ 


tation). 

In conclusion, with our study of 45 member stars in 
NGC 6139 conducted with FLAMES spectra we determined 
for the first time its metallicity with high-resolution spec¬ 
troscopy (the cluster has an intermediate metallicity, with 
[Fe/H]=-1.58 dex). We also measured the light elements O, Na, 
Mg, and Al involved in proton-capture reactions of H-burning 
at high temperatures; we detected the usual correlations and an¬ 
ticorrelations found in almost all GCs investigated so far. The 
ratio of hrst-to-second generation stars in NGC 6139 is typi¬ 
cal of the majority of GCs and is consistent with its relatively 
high total present-day mass. We do not support the idea that 
NGC 6139 is a FG-dominated cluster. The intermediate exten¬ 
sion of the Na-O anticorrelation, as measured by the interquar¬ 
tile range of the [O/Na] ratio, fits very well the relation with the 


maxim um temperature along the HB found by ICarretta et al 


(2007)- All in all NGC 6139 behaves like a normal MW GC 
of similar mass and HB extension. The formation event that 
shaped its internal chemistry seems consistent with the pattern 
of a regular, relatively high-mass, blue HB GC. 

However, the finding of a large variation in Al, which 
should also imply some spread in He, does not perfectly fit 
with the HB of the cluster. We calculated IQR[A1/Mg], finding 
a value of 0.425; according to Fig. 28 in Gratton et al. ( 2Q10l) . 
this would require a AY of about 0.02 (but note that the rela¬ 
tion between IQR[A1/Mg] and AY is rather poorly determined. 


being based on only eight points). We note that also the more 
massive, meta l-poore r NGC 5024 (M53, with My - -8.71, 
[F e/H]=-2.10, Harrisl 19961). proposed as a FG-dominated GC 
by Caloi & D Antonal d201 ll) on the basis of its HB and dis¬ 
cussed at length in tha t pap er, lias been demonstrated to host a 
substantial SG by Meszaros et al.l ( 2015 ), using APOGEE data. 
NGC 5024 has a large spread in Al (with a bimodal distribu¬ 
tion) and almost no Mg variation; it_behaves like NGC 5272 
(M3, M v = -8.88, rFe/Hl=-1.50. lHarrislll996l) . another GC 
with a rather short HB, even more massive than NGC 6139, 
but with a similar metallicity. The relation between Al and He 
is evidently not com pletely understood. As as an additionally 


cautionary note, Bastian et al.l (2015) claim that no one of the 


enrichment mechanisms proposed to explain multiple popula¬ 
tions in GCs can consistently do so, because all predict too 
large He spreads in order to reproduce the observed spreads 
in light elements, notably those in Na and O. More observa¬ 
tional constraints are required, both from spectroscopy with 
large samples of stars for which elements of all nucleosynthetic 
chains are measured (like we and other groups are doing for 
an increasing number of of GCs) and from precise photometry 
(e.g., Piotto et al. . 20151 to determine He). These constraints 
will be helpful to improve modelling the formation and early 
enrichment phases of GCs. 
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Table 2. Information on the member stars observed, the complete table is available in the on-line version 


ID 

RA 

(hh:mm:ss) 

Dec 

(dd:pp:ss) 

V 

B 

K 

(2MASS) 

RV 
(km s -1 ) 

err 

(km s" 1 ) 

Notes 

t00540 

16:27:39.82 

-38:50:18.05 

15.825 

UVES 

17.656 

10.601 

24.00 

0.30 

4, RVji2=37 

100553 

16:27:39.22 

-38:51:18.04 

15.859 

17.557 

10.856 

29.63 

0.60 

1, RV S12 =46 

t00670 

16:27:42.77 

-38:50:56.29 

16.069 

17.703 

11.329 

21.27 

0.50 

3, RVsi2=36, AGB 

100951 

16:27:35.86 

-38:53:40.28 

16.391 

18.077 

11.706 

26.21 

0.60 

13, RVs 12=44 

tO 1348 

16:27:40.82 

-38:48:42.10 

16.764 

18.379 

12.049 

35.00 

0.30 

11, RV S12 =19 

101554 

16:27:38.41 

-38:49:23.44 

16.920 

18.434 

12.261 

26.71 

1.60 

8, RV S12 =38, AGB 

101699 

16:27:42.55 

-38:49:53.35 

17.010 

18.605 

12.383 

26.09 

0.50 

7,RV S12 =42 

t00261 

16:27:43.32 

-38:57:04.20 

GIRAFFE, members 
15.097 17.217 9.530 

17.37 

0.35 


t00284 

16:27:46.06 

-38:54:03.74 

15.218 

17.170 

9.988 

21.40 

0.55 


t00316 

16:27:46.75 

-38:51:04.61 

15.337 

17.240 

10.145 

15.68 

0.34 


t00329 

16:27:20.70 

-38:53:22.86 

15.341 

17.445 

9.744 

28.93 

0.42 


t00337 

16:27:30.42 

-38:52:04.70 

15.363 

17.450 

9.754 

25.63 

0.84 



In Notes we give the number and RV in Saviane et al. (2012); 
Star t06169 is on the HB. 


Table 3. Information on non member stars on the basis of RV or metallicity, all observed with GIRAFFE, the complete table is 
available in the on-line version 


ID 

RA 

(hh:mm:ss) 

Dec 

(dd:pp:ss) 

V 

B 

K 

(2MASS) 

RV 
(km s -1 ) 

err 

(km s" 1 ) 

Note 

tOO187 

16:27:25.45 

-38:48:42.66 

14.781 

17.020 

7.078 

-16.98 

1.00 

RV 

t00306 

16:27:10.97 

-38:53:19.81 

15.271 

17.381 

9.388 

-85.72 

0.49 

RV 

t00392 

16:27:48.91 

-38:51:43.38 

15.555 

17.294 

10.639 

-34.11 

0.99 

RV 

t00466 

16:28:07.92 

-38:49:47.99 

15.666 

17.730 

9.262 

-181.62 

0.69 

RV 

100513 

16:28:00.26 

-38:52:21.00 

15.737 

17.832 

9.739 

-116.30 

0.48 

RV 

100526 

16:27:45.67 

-38:55:14.59 

15.764 

17.815 

9.802 

-160.69 

0.52 

RV 


Table 4. Adopted atmospheric parameters and derived metallicity for confirmed member stars. the complete table is available 
in the on-line version 


Star 

7eff 

(K) 

l°gg 

(dex) 

[A/H] 

(dex) 

v, 

(km s~ 

nr 

*) 

[Fe/H]i 

(dex) 

rms 

nr 

[Fe/Hn 

(dex) 

rms 

100540 

4156 

0.95 

-1.55 

1.85 

UVES 

121 

-1.545 

0.089 

20 

-1.551 

0.076 

100553 

4212 

1.07 

-1.61 

1.40 

107 

-1.610 

0.102 

20 

-1.505 

0.131 

100670 

4417 

1.28 

-1.57 

1.43 

91 

-1.570 

0.112 

13 

-1.566 

0.164 

t00951 

4398 

1.44 

-1.61 

1.93 

83 

-1.614 

0.095 

17 

-1.475 

0.098 

101348 

4473 

1.57 

-1.57 

1.52 

74 

-1.571 

0.101 

15 

-1.586 

0.109 

101554 

4634 

1.76 

-1.63 

0.98 

52 

-1.629 

0.140 

10 

-1.588 

0.173 

101699 

4546 

1.72 

-1.52 

1.64 

76 

-1.517 

0.090 

15 

-1.517 

0.105 

t00261 

3922 

0.51 

-1.56 

1 

2.09 

GIRAFFE 

40 

-1.558 

0.087 

3 

-1.524 

0.207 

t00284 

4022 

0.71 

-1.59 

1.97 

36 

-1.592 

0.099 

3 

-1.517 

0.064 

100316 

4056 

0.77 

-1.58 

2.03 

48 

-1.576 

0.123 

3 

-1.576 

0.216 

100329 

3968 

0.59 

-1.59 

2.02 

38 

-1.588 

0.079 

3 

-1.514 

0.089 

100337 

3971 

0.59 

-1.61 

2.08 

41 

-1.605 

0.110 

3 

-1.512 

0.016 
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Table 5. Light element abundances for confirmed member stars, the complete table is available in the on-line version 


Star 

nr 

[0/Fe]i 

rms 

nr 

[Na/Fe]i 

rms 

nr 

[Mg/Fe]i 

rms 

nr 

[Al/Fe]i 

rms Lim? 

PIE 







UVES 







t00540 

i 

-0.041 


3 

0.637 

0.010 

3 

0.518 

0.068 

2 

0.802 

0.008 

I 

t00553 

2 

0.174 

0.071 

3 

0.303 

0.050 

3 

0.459 

0.038 

2 

0.395 

0.084 

I 

t00670 

1 

-0.126 


2 

0.377 

0.089 

2 

0.420 

0.071 

2 

0.559 

0.006 

I 

t00951 

1 

0.506 


1 

-0.028 


2 

0.482 

0.011 

2 

0.330 

0.126 

P 

t01348 

1 

0.369 


2 

0.030 

0.055 

2 

0.464 

0.006 

1 

0.270 


P 

101554 

1 

0.263 


2 

0.313 

0.061 

2 

0.532 

0.204 

2 

0.724 

0.097 

I 

101699 

1 

0.146 


2 

0.414 

0.001 

2 

0.514 

0.102 

2 

0.536 

0.112 

I 







GIRAFFE 







t00261 

2 

0.078 

0.185 

4 

0.631 

0.040 

3 

0.503 

0.072 

2 

0.555 

0.031 

I 

t00284 

2 

0.209 

0.006 

4 

0.657 

0.040 

3 

0.563 

0.057 

2 

0.556 

0.000 

I 

t00316 

1 

0.037 


4 

0.637 

0.061 

2 

0.503 

0.006 

2 

0.940 

0.128 

I 

t00329 

2 

0.430 

0.051 

4 

0.340 

0.036 

3 

0.538 

0.015 

2 

0.279 

0.006 

I 

100337 

2 

0.105 

0.006 

4 

0.747 

0.019 

3 

0.516 

0.076 

2 

0.921 

0.040 

I 





